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Abstract 

We present a microscopic model for the calculation of a-decay half lives employing 
potentials obtained from relativistic and non-relativistic self-consistent mean-field 
models. The nuclear and Coulomb potentials are used to obtain the tunneling prob¬ 
ability and, in one model variant, also the knocking frequency. The model contains 
only one parameter. We compare this approach employing several modern mean- 
field parametrizations to experimental data and to the semi-empirical Viola-Seaborg 
systematics. We extrapolate our model to superheavy nuclei where assumptions en¬ 
tering semi-empirical approaches might lose validity. 
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1 Introduction 


Alpha decay was experimentally studied by Rutherford just before the begin¬ 
ning of the 20th century and constitutes a beautiful example of quantum me¬ 
chanical tunneling through a (classically impenetrable) barrier. The o-particle 
or helium nucleus, which is formed within the mother nucleus, tunnels through 
the Coulomb barrier created by the protons. In 1928, after the advent of quan¬ 
tum mechanics, Gamow proposed a semi-classical interpretation for o-decay 
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[1], His model implied that the a-decay constant A can be written as a prod¬ 
uct of three factors: the probability / that an alpha particle is formed inside 
the nucleus {pre-formation factor), the frequency z/ at which the alpha parti¬ 
cle knocks against the nucleus’s potential well {knocking frequency), and the 
probability P with which the particle is able to tunnel through the barrier 
{Gamov factor), leaving behind the {Z — 2, N — 2) daughter nucleus. Thus, in 
this model the a decay constant is given by 

\ = f xux P (1) 

We note that this approach assumes the statistical independence of these three 
processes. 

The probability of transmission through the barrier, P = e’’’, can be deter¬ 
mined by calculating semi-classically in the WKB approximation [2] 



The knocking frequency u is given by [3] 



which, for a constant nuclear potential V, reduces to v/2R, i.e., the number of 
times a particle with velocity v (determined from its energy E) goes back and 
forth in a conhned space of size 2R. In these formulas, p is the reduced mass 
of the daughter-alpha particle system, and ri and r 2 are the points where the 
energy of the alpha particle is equal to the potential, ri is the inner, r 2 is the 
outer turning point, respectively. 

Eq. (3) takes into account the varying contribution of the kinetic to the total 
a-particle energy which depends on the potential depth at position r. Time- 
dependent quantum-mechanical calculations [3], however, indicate that the 
oscillations of the a-particle wave-function within the nucleus have an ex¬ 
tremely small amplitude, which might lead to a better estimate of the number 
of barrier assaults per second. 


The a-decay half live Ti /2 is then hnally obtained by 

In 2 

r ./2 - — 


( 4 ) 


A semi-empirical formula for the calculation of a half lives was presented by 
V. E. Viola and G. T. Seaborg (VS) in 1966 [2]. It is based on the Gamow 
picture and assumes a nuclear potential of uniform depth that rises abruptly 
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at the radius of the nucleus, going from this uniform potential to a Coulomb 
potential, resulting in a cusp at the radius of the nucleus. Also, the nucleus 
and therefore the potential are assumed to be spherically symmetric. 

The resulting formula reads 


logri/2 


aZ + b 

y/Qa 


+ (^cZ + d), 


( 5 ) 


yielding the a-decay half life Ti /2 in units of seconds. The experimental input 
are the charge number Z and the Qa value (in units of MeV) for a given nucleus 
{Z, N). It contains four adjustable parameters a, b, c, d. These parameters have 
been adjusted to the experimental data on half lives known at that time. In 
more recent work, the parameters have been readjusted, taking into account 
new experimental data on superheavy elements [4]. Table 1 summarizes the 
values of these parameters and also contains the parameters of our adjustment 
made in this paper, see Subsection 3.1 for details. The VS systematics has 
been used in Ref. [5] to calculate a half lives for unknown superheavy nuclei. 
We note that there exists a variety of very accurate semi-empirical half live 
formulas, see for example Refs. [6,7]. 



a 

b 

c 

d 

original VS 

2.11329 

-48.9879 

-0.390040 

-16.9543 

readjusted VS 

1.66175 

-8.5166 

-0.20228 

-33.9069 

our adjustment 

1.58134 

-2.27889 

-0.23543 

-30.1503 


Table 1 

Parameters of the original VS systematics [2], its readjustment including data on 
superheavy nuclei [4], as well as our adjustment to experimental data from Ref. [8]. 


In this paper, we present a model for the calculation of a-decay half lives based 
on microscopically calculated potentials from self-consistent mean-held mod¬ 
els, namely the relativistic mean-held (RMF) model and the Skyrme-Hartree- 
Fock (SHF) approach, thus replacing the simplihed assumptions of Eq. (5). 
In Section 2 we present our approach. Section 3 describes the adjustment and 
analysis of our method. We compare its predictions with both experimental 
data as well as a new adjustment of the VS systematics. In Section 4 we ex¬ 
trapolate to superheavy nuclei and compare the two approaches. We conclude 
in Section 5. 
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2 Theoretical Framework 


In our approach, which we henceforth denote as o-mf, the o-decay constant 
is given by Eq. (1). The o-nucleus potential is given by 

Va = 2xVp + 2xVn ( 6 ) 

The proton potential Vp and the neutron potential Vn are taken from spherical 
self-consistent mean-held calculations (see below). Thus, Va contains the sum 
of the nuclear and the Coulomb potential. In the relativistic models, the Va 
potential is obtained as the sum of the scalar and vector potentials. The a- 
particle is treated as a spin-saturated system of two protons and two neutrons, 
thus it does not feel the spin-orbit force that is generated by the strong scalar 
and vector helds adding up with the same sign in the relativistic approaches, 
and put in by hand as an additional term in the SHF approach. 

The transmission probability and the knocking frequency, given by Eqs. (2) 
and (3), respectively, are calculated by numerical integration using the mean- 
held potentials. While the bosonic and structure-less a-particle picture can be 
expected to work for the tunneling process as long as the Compton wavelength 
of the a-particle is small compared to the width of the barrier, the ’knock¬ 
ing’ process of the a-particle within the Coulomb barrier is certainly more 
complicated due to the fermionic nature of the nucleons and the necessary 
anti-symmetrization between the wave functions of a-particle and nucleus. 
Still, we may expect that these corrections will not ahect too much the overall 
performance of our model. 

Our approach contains only one free parameter which we denote by / and 
which is related to / in Eq. (1). It needs to be htted to experimental data. 
This parameter not only in an average way describes the formation probability, 
but also absorbs all errors and approximations made in the calculations of P 
and u. Interpreting it as an (average) probability, it should be less than or 
equal to 1. According to Ref. [9] a pre-formation factor should be below 0.1. 
In order to estimate the accuracy of the knocking frequency calculation, we 
also build combined parametrizations with the combined parameter c = f x u, 
in which only the transmission probability is computed numerically. 

We note that with the knowledge of the single-particle wave functions, it is 
possible to calculate the pre-formation factor, see for example Ref. [10] for a 
calculation using Nilsson model single particle wave functions and Ref. [11] 
for a parameter-free approach for calculating the a-decay half live of ^^^Po 
employing a combined shell and cluster model. Complementing our model by 
such a pre-formation calculation would lead to a fully microscopic, parameter- 
free approach, which is, however, beyond the scope of this paper and will 
be dealt with in forthcoming work. The assumption of spherical potentials 
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in our approach loses its validity for deformed nuclei. It is, however, also 
used in the Viola-Seaborg systematics and appears to work well for life-time 
predictions. Future work will include potentials taken from deformed mean- 
held calculations [12]. 

The mean-held a-potential corrects the simplihed potential of other approaches 
by providing a smooth surface and r-dependent interior nuclear potential, and 
eliminates the spike at the transition from the nuclear to the Coulomb part, 
see Fig. 1 (the authors of Ref. [13] use a non-self-consistent potential with 
a smooth surface). Also, its radius automatically possesses the correct mass- 
dependence. The accuracy of our model tests the asymptotics of the Coulomb 
contribution to the mean-held potentials stemming from the relativistic mean- 
held calculations up to typical radii of r ^ 70 fm. The mean-held potentials 

- due to their self-consistent interdependence of proton and neutron densities 

- automatically introduce a dependence on neutron number, i.e., an isospin 
dependence. For each nucleus, the input to our model is the Vk potential from 
the daughter nucleus (assuming that the a particle has already formed), the 
Qd value and the mass of the nucleus. The potential is being represented on 
a radial 1-D grid in coordinate space. The classical turning points ri and 
are determined and then used as boundaries for the numerical computation of 
transmission and knocking frequency according to Eqs. (2) and (3). Because 
of the uncertainties related to the calculation of odd-even and odd-odd nuclei, 
we consider only even-even nuclei in this study. 

The self-consistent mean-held models used in this work are on the one hand 
two variants of the relativistic mean-held model, one employing contact in¬ 
teractions (RMF-PC) , and the other variant employing hnite-range boson 
helds (RMF-FR), and on the other hand the Skyrme-Hartree-Fock (SHF) ap¬ 
proach. We refer the reader to Refs. [14,15,16,17] for a detailed discussion of 
these models. In modern terms, self-consistent mean-held models are approx¬ 
imations to the exact many-body functional in the spirit of density functional 
theory [18], their interaction terms governed by principles of ehective held 
theory such as naive dimensional analysis. Thus, these models are able to ab¬ 
sorb various ground-state correlations and to incorporate physics beyond the 
(literal) mean-held approximation [19]. 

The models used here contain between 6-12 parameters plus two pairing 
strengths (for protons and neutrons separately) for BCS pairing with a density- 
independent (5-force. The parameter sets employed in this work are NL3 [20] 
and NL-Z2 [21] for RMF-FR, PC-Fl [22] for RMF-PC, and SkI3 [23] and 
SLy6 [24] for SHF, which are modern mean-held forces that deliver great pre¬ 
dictive power. They result from careful adjustments to nuclear ground-state 
observables, e.g., binding energies, rms charge and dihraction radii, and sur¬ 
face thicknesses [22]. The nucleon single-particle wave-functions are calculated 
self-consistently on a grid in coordinate space employing matrix multiplica- 
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tions in Fourier space for the derivatives. To obtain the solution with minimal 
energy, the damped gradient-step method is used [25]. Calculations with both 
models share the same basic numerical routines, thus numerical differences 
do not interfere with differences in physics for model comparisons. The center 
of mass correction is used as in the adjustment procedure of the mean-held 
parametrization. 

Since the a-particle has a hnite extension, in addition we test another model 
variant by folding its parametrized density distribution with the mean-held 
alpha potential, i.e. 


= I p^{f- R) f4(f) d^r (7) 

The radial density parametrization of the a-particle is taken from Ref. [26] 
and is given by 

p„(r) = 0.4229 X (8) 

where Pa{r) is given in units of fm“^, and the radius r needs to be specihed 
in units of fm. The ehect of the folding procedure on the mean-held potential 
is shown in Figure 1. 

Our a-decay model contains no assumptions on the nuclear potential, except 
for the calculation of the potentials in spherical symmetry. Thus, provided 
that the potentials from the self-consistent mean-held models are realistic, it 
should lead to more reliable extrapolations to exotic and superheavy nuclei 
compared to semi-empirical models. We note that these potentials automat¬ 
ically possess - due to the adjustment of the mean-held models to density- 
related ground-state observables - the correct mass dependence of the poten¬ 
tial radius and also information on the surface region of the nucleus. These 
need to be parametrized in more schematic approaches. 

We would like to point out possible extensions of the a-mf approach, most 
of them going clearly beyond the scope of semi-empirical approaches. The 
microscopic calculation of the pre-formation factor has already been men¬ 
tioned. Another enhancement would be the usage of 2-dimensional potentials 
from axially deformed nuclei. In the case of shape coexistence of nuclei, the 
dependence of the decay width on the prolate or oblate minimum can be in¬ 
vestigated. Furthermore, decay widths stemming from the ground state or an 
isomeric state (shape isomer) of a heavy nucleus can be computed. The de¬ 
scription of hne structure in a-decay will be an interesting application, see 
for example Ref. [27]. The model could be enhanced still further by replac¬ 
ing the WKB approximation with numerical calculations of the stationary or 
time-dependent Schrbdinger equation. 
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Fig. 1. Alpha-nucleus potential obtained from a mean-field calculation with PC-Fl 
(dashed line), the same potential but folded with the a-particle density (dotted 
line), and a schematic potential according to the VS systematics (full line) for the 
nucleus ^s^cf 

3 Results 


3.1 Adjustment 


For the model adjustment and comparison, we use recent a-decay data from 
Ref. [8]. Since the experimental data varies over thirty orders of magnitude 
(from 1.1 X 10“^ s to 2.2 x 10^^ s) we have performed a logarithmic fit to 
adjust the model parameter. The logarithmic error of a given data point is 
given by 


X = 



Tll2,th 

R/2,ea;p 


( 9 ) 


Computing errors in this way measures how many orders of magnitude the 
calculated point is away from the experimental value (for example, y = 2.5 
means that the calculated value is 2.5 orders greater than the experimental 
value, while y = —3 means that the calculated value is 3 orders of magnitude 
less than the experimental value). 


To obtain a fair comparison between a-mf and VS systematics, we have read¬ 
justed the VS systematics to exactly the same data that has been used in the 
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adjustment of our model [8], employing the Monte-Carlo methods described in 
Ref. [22] (down-hill methods got stuck very easily in adjacent local minima). 
The Monte-Carlo optimization was finally followed by a down-hill minimiza¬ 
tion, driving the parameter vector to the bottom of the (local) minimum. The 
resulting parameters are shown in Table 1. Note that this adjustment has 
been started with the parameters of the readjusted VS formula. The usage 
of Monte-Carlo techniques has reduced considerably. However, there might 
very well exist several local minima, of which some could correspond to models 
with even better predictive power. 


The values of the a-mf approach employing different mean-field forces are 
displayed in table 2. In the following we implicitly refer to x^ per point when 
we write x^■ The x^ value of our VS systematics adjustment is included as well. 
We see that the Skyrme forces SkI3 and SLy6 in three of four cases deliver a 
X^ value that is below x^ of this adjustment of the VS systematics, while the 
relativistic models have slightly larger values. This model difference is quite 
interesting since it yields information on the (otherwise unobservable) poten¬ 
tials of the different approaches. The SHF potentials (for both Skyrme forces) 
lead to better predictive power compared to RMF. This model difference, that 
might be related to the effective masses, radial dependences, and Coulomb ex¬ 
change (absent in RMF) deserves further attention in future investigations. 
The model variants with the combined parameter c (absorbing both prefor¬ 
mation factor and knocking frequency into one parameter) perform better in 
most cases. The differences are not dramatic, however, indicating that the ex¬ 
plicit calculation of v is consistent with our model assumptions. The (average) 
preformation factors / are smaller than one and thus can be interpreted as 
probabilities. They also agree with the estimates from Ref. [9]. 


The model variant employing the folded potentials from calculations with 
PC-Fl performs slightly better when explicitely computing the knocking fre¬ 
quency, and slightly worse when absorbing it into the parameter c. This might 
be connected to the fact that in the folding case we still use the simple reduced 
mass in the calculation of the tunneling probability, while the mass parameter 
is certainly different for an extended object. The parameter / for the folded 
PC-Fl potential is a factor of 10 smaller than the value without folding, which 
is related to the increase of the tunneling probability in the folded case, see 
Fig. 1. 


We have determined the x^ values of the previous hts of the VS systematics 
using the experimental data from Ref. [8]. The original Viola-Seaborg method 
yields x^ = 2.57, and the improved Viola-Seaborg method x^ = 1.27. 



force 

/ 


c (1/s) 


NL3 

0.028 

0.21 

9.78el9 

0.21 

NL-Z2 

0.014 

0.17 

4.57el9 

0.16 

PC-Fl 

0.080 

0.24 

3.09e20 

0.22 

PC-Fl, folded 

0.006 

0.23 

1.95el9 

0.24 

SkI3 

0.087 

0.17 

3.46e20 

0.15 

SLy6 

0.075 

0.15 

2.85e20 

0.14 

our VS fit 


0.16 




Table 2 

The parameter / and corresponding values for the model variant with explicit 
calculation of the knocking frequency (columns 1 and 2), and the combined pa¬ 
rameter c and corresponding values for the combined variant (columns 3 and 
4), for the mean-field forces as indicated. For comparison, the x^ value of our VS 
systematics adjustment is shown. 

3.2 Analysis 


In Figure 2 we compare half lives between a-mf employing the relativistic 
mean-held parametrization PC-Fl and experimental data (since the differ¬ 
ent mean-held approaches have similar predictive power, comparisons of other 
forces with experiment yield similar hgures). Note that the a-mf approach 
delivers a very good overall agreement with |x| < 0.3 — 0.6 in most of the 
cases, i.e., the deviations are well below an order of magnitude. These errors 
have predominantly positive signs corresponding to an over-estimation of the 
a life-times, i.e., an over-estimation of stability of these nuclei with respect to 
a-decay. Negative errors are visible at and near shell closures, see the discus¬ 
sion below. The accuracy of our approach does not appear to depend strongly 
on the mass number of the decaying nucleus, indicating that the mass num¬ 
ber dependence is correctly reproduced in the mean-held potentials. The best 
agreement, however, is achieved for the actinides. 

We would like to examine our predictions for indications of physics not (yet) 
accounted for in our model. We can expect to see deviations originating from 
i) shell ehects and ii) deformation ehects. Let us hrst examine shell ehects. At 
the magic neutron number N = 126, for Z = 84 —90 isotopes we notify a jump 
in the logarithmic error: a-mf overestimates the stability of = 128 isotones, 
while the stability of = 126 isotones is underestimated (indicated by a 
negative error). We see a similar ehect at the magic proton number Z = 82: a- 
mf underestimates stability. This trend can be nicely observed for the N = 128 
isotones. As the proton number drops oh towards 82, the nuclides become 
more and more stable compared to our predictions. Even though shell ehects 
are taken into account to some extent via the experimental Qa values, the 
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Fig. 2. Logarithmic error of our approach with PC-Fl as a contour plot versus 
neutron number (x-axis) and proton number (y-axis). Each square corresponds to 
an even-even nucleus. 

computation of the tunneling probability cannot reproduce them sufficiently. 
They would enter, however, in a microscopic calculation of the preformation 
amplitude. 

Another cause of trends is the deformation of nuclei. Our current model - 
as does the VS systematics - assumes spherically symmetric nuclei. It is in¬ 
teresting that the assumption of spherical potentials works so well both in 
semi-empirical approaches and in our model. However, most nuclei for which 
we predict a half lives are deformed. Thus, for axially deformed nuclei the 
Coulomb barrier is 2-dimensional, leading to a more complicated tunneling 
dynamics (see Refs. [28,29] for a time-dependent description of 2-dimensional 
tunneling in proton emission). For the most part, our predictions do not dis¬ 
play a random distribution pattern, but rather show clusters where the model 
is more accurate and clusters were it is less accurate. For example, nuclei with 
atomic numbers close to Z = 88 and N = 136 are less stable than our model 
predicts, and around Z = 102 and N = 152 there is a region of extra stability 
that our model does not yet account for. It is not possible to tell directly from 
our data whether the deformation effect results in making nuclei appear to 
have extra stability or extra instability, due to the masking effect of optimiz¬ 
ing the (average) pre-formation factor. We believe, however, that accounting 
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Fig. 3. Logarithmic error of the newly adjusted Viola-Seaborg systematics as a 
contour plot versus neutron number (x-axis) and proton number (y-axis). Each 
square corresponds to an even-even nucleus. 


for deformations in nuclei shall improve our model. 

Comparing the predictions of the a-mf approach (Fig. 2) with the ones from 
our VS systematics adjustment (Fig. 3), we detect a few differences. Firstly, 
nuclei with larger errors with respect to the surrounding ones are distributed 
differently in both approaches. Secondly, there are small regions of nuclei which 
are described better in each model, for example the nuclei centered around 
Z = 86, iV = 130 have larger errors within the VS systematics, while nuclei at 
or next to the a-decay chain starting aX Z = 90, N = 120 are more accurately 
described by the VS model. Moreover, we detect similar wrong trends that 
can be attributed to missing deformation and/or shell structure features in 
the VS systematics. 

In the next Section, we extrapolate with a-mf employing SLy6 to superheavy 
nuclei and compare with the predictions of the VS systematics. 
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Fig. 4. Logarithmic ratio of the predicted life time from a-mf over the predicted 
life time of our Viola-Seaborg adjustment systematics as a contour plot versus neu¬ 
tron number (x-axis) and proton number (y-axis). Each square corresponds to an 
even-even nucleus, nuclei without a converged solution have been omitted. The black 
squares indicate the valley of (3 stability, the grey line corresponds to the two-proton 
dripline. 

4 Extrapolation to superheavy elements 


We employ the a-mf approach with the force SLy6, which has performed 
best in terms of for an extrapolation to superheavy nuclei. We compare its 
predictions to the predictions from our adjustment of the VS systematics. The 
Qa values are obtained from axially deformed ground-state calculations with 
SLy6, similarly to Refs. [30,5]. The potentials are obtained from calculations 
in spherical symmetry. For both models, the Qa values calculated with SLy6 
are used, since a) experimental data in that region is very scarce (especially 
for even-even nuclei) and b) we aim at a pure model comparison. 

The logarithmic ratio of the two approaches can be seen in Fig. 4. The differ¬ 
ence between the two approaches remains quite small (below 0.1) for nuclei 
close to /3-stability. Also, both models show no visible difference in the error 
with mass number, a feature we have seen before for the known elements. 
However, an isovector trend appears in the hgure between the two methods. 
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With decreasing N — Z ratio, the a-mf approach predicts increasingly larger 
life times. Since both approaches use the same - calculated - Qa values, the 
isospin dependence has its origin in the potential shapes used in a-mf which 
introduces a dependence on the neutron number. This isovector trend is absent 
in the VS systematics. 

For nuclei with masses beyond superheavy nuclei {A > 292), various groups 
have theoretically found semi-bubble and bubble structures [31,32,33,34], i.e., 
the density in the interior of the nucleus is reduced or even completely sup¬ 
pressed, leading to nuclear systems with an inner surface. Because of this 
structure, these systems possess potentials that differ to a great extent from 
usual nuclear potentials. Thus, they will constitute an interesting application 
for a-decay models based on self-consistent potentials [12]. 


5 Conclusions 


We have presented a one-parameter model for o-decay that can be used to 
predict life times of o-emitters. This model can be applied successfully over a 
large range of nuclei. Its accuracy is comparable to or even slightly better than 
semi-empirical approaches. This has been demonstrated by a new adjustment 
of the Viola-Seaborg formula. However, universality and extensibility are the 
decisive features of this model. It can be applied to exotic and superheavy 
nuclei, in which the assumptions of (very successful) semi-empirical formulas 
about the shape of the potential lose more and more validity, given that the 
potentials from the self-consistent mean-held models are still realistic. By em¬ 
ploying a variety of recent mean-held parametrizations, we have found that 
a-nucleus potentials from the Skyrme-Hartree-Fock method yield slightly bet¬ 
ter results compared to the RMF model potentials. 

An extrapolation to superheavy elements has shown that the a-mf approach 
displays an isospin dependence of life times which is absent in the Viola- 
Seaborg approach. This point deserves further attention as new data become 
available. 

The model is simple in terms of required parameters. While the Viola-Seaborg 
systematics contains four parameters that need to be simultaneously opti¬ 
mized, it contains only one adjustable parameter, the (average) pre-formation 
(or a combined) factor. Furthermore, our model has ample room for im¬ 
provement. One potential improvement would be to microscopically calculate 
the pre-formation factors using nuclear single-particle and a-particle wave- 
functions, yielding a fully microscopic, parameter-free approach. Other en¬ 
hancements and applications involve multi-dimensional potentials and hue- 
structure a-decay calculations. 
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